Non-Boltzmann Modeling for Air Shock-Layer Radiation at Lunar- Return Conditions 
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Abstract 

This paper investigates the non-Boltzmann modeling of the radiating atomic and molecular 
electronic states present in lunar-return shock-layers. The Master Equation is derived for a 
general atom or molecule while accounting for a variety of excitation and de-excitation 
mechanisms. A new set of electronic-impact excitation rates is compiled for N, O, and N 2 + , 
which are the main radiating species for most lunar-return shock-layers. Based on these new 
rates, a novel approach of curve-fitting the non-Boltzmann populations of the radiating atomic 
and molecular states is developed. This new approach provides a simple and accurate method for 
calculating the atomic and molecular non-Boltzmann populations while avoiding the matrix 
inversion procedure required for the detailed solution of the Master Equation. The radiative flux 
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values predicted by the present detailed non-Boltzmann model and the approximate curve-fitting 


approach are shown to agree within 5% for the Fire 1634 s case. 
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Nomenclature 

Bohr radius, equal to 0.529x10 8 cm 

transition probability for a line with an upper level j ’ and a lower level i ’ 

velocity of light, equal to 2.997925xl0 10 cm/s 

electronic term energy for an atomic level i (cm’ 1 ) 

ionization energy of an atomic level i 

ionization energy of hydrogen, equal to 109,697 cm" 1 

electron charge, equal to 4.80298xl0" 10 cm 3/2 g 1/2 /s 

degeneracy for an atomic level i 

Planck’s constant, equal to 6.6256x10 27 erg-s 

frequency integrated emission coefficient from atomic lines (erg/cm /s/sr) 

Boltzmann constant, equal to 1.38054x10 16 erg/K 

electron impact excitation rate from level i to level j; the presence of a d or c in 
parenthesis refers to a dissociated or ionized level, respectively, 
heavy-particle impact excitation rate from level i to level j; the presence of ad or c 
in parenthesis refers to a dissociated or ionized level, respectively. 

no 

electron mass, equal to 9.1091x10' g and the number of electronic states treated for 
a given atom. 

molecular weight of species i (kg/kg-mol) 
number density of an atom 

number density of an atomic level i (particles/cm ) 
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N+ = ion number density (particles/cm ) 

N e = electron number density (particles/cm ) 

q, - = wall-directed radiative heat flux, calculated in units of erg/s/cm 2 , but usually 

presented in units of W/cm 2 (1x10 7 W/cm 2 = lerg/s/cm 2 ) 

Q x = partition function of the energy mode x or an atom a (nondimensional) 

r 0 = vector defined in Eq. (13) for the Saha decrements in Eq. (12) 

r 00 = vector defined in Eq. (13) for the Saha decrements in Eq. (12) 

r . = vector defined in Eq. (13) for the Saha decrements in Eq. (12), 

where j represents a level less than or equal to the number of levels that the QSS 
assumption is applied to. 

T e = electronic temperature (K), assumed equal to T v and T e 

x = number of active electrons for an atomic level 

Ajj = escape factor for the transition from i to j (nondimensional) 

Pi = refers to the Saha-decrement of level i, defined by Eq. (8) 

t, = radiative lifetime for a transition with an upper state i (s) 

Subscripts 

e = refers to the electronic energy mode 

EQ = assumes that Boltzmann and Saha equilibrium exists among the electronic levels 
i = refers to the lower electronic state of an atomic transition 

j = refers to the upper electronic state of an atomic transition 

Superscripts 

SB = indicates Saha-Boltzmann equilibrium 

= indicates the radiative flux or intensity directed towards the vehicle wall 
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I. Introduction 


In a companion paper 1 the calculation of the radiative emission and absorption is discussed for 
the various radiative processes of atoms and molecules. The magnitude of the emission and 
absorption coefficients from these various radiative processes are proportional to the number 
density of the emitting or absorbing electronic state. For all of the examples worked by Johnston 
et al. 1 , it was assumed that these electronic states are populated by a Boltzmann distribution. This 
assumption is accurate for chemical equilibrium conditions. It has been shown' that there are 
significant regions of chemical nonequilibrium in the lunar return shock-layers of present 
interest, which require that the non-Boltzmann electronic state populations be considered. The 
present paper investigates the calculation of these non-Boltzmann electronic state populations for 
the significant radiating species present in air shock-layers at lunar return conditions. The species 
considered in this paper include two atomic species, N and O, and one molecular species, N? + . 
Although it was shown by Johnston et al . 1 that the NO and O2 molecules contribute to the 
radiation at lunar-return conditions, it was found that they could be assumed to be populated in a 
Boltzmann distribution. This assumption is allowable because these molecules do not reach a 
peak in the nonequilibrium region of the shock layer (unlike the N^ 4 * molecule), so their 
nonequilibrium contribution is small. 

A review of past non-Boltzmann modeling for shock-layer radiative heating is discussed 
in Section II. The definition of the atomic levels for N and O used in the present work are 
discussed in Section III. The atomic and molecular electronic state-populations are modeled in 
the present work using a collisional-radiative (CR) model, which accounts for the repopulation 
and depopulation of the atomic and molecular states through collisional and radiative processes. 
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These processes are discussed individually in Section IV. In Section V, the general formulation 
of the CR model’s governing equation, called the Master Equation, is derived from the processes 
discussed in Section IV. A practical method of solving the Master Equation is presented in 
Section VI. The task of assembling the rates required for the atomic CR modeling of N and O is 
discussed in Section V. Excitation rates from several sources are compared and a final set of 
rates is compiled and presented. Based on this detailed rate model and the solution procedure for 
the Master Equation discussed in Section VI, an “approximate atomic CR” (AACR) model is 
developed, which allows for the efficient and accurate application of the atomic CR model. The 
AACR model is presented in Section VIII, and is conveniently represented by curve fits 
presented in Appendix C of Ref. 18. The non-Boltzmann modeling of the N 2 + is discussed in 
Section IX. The influence of various collisional processes on this molecule is studied, and rates 
from several sources are compared. A recommended set of rates is proposed, and a set of curve 
fits are constructed for the radiating electronic states of N 2 + . 

II. Overview of Non-Boltzmann Modeling 

A significant difference between two widely used radiation codes is the treatment of the 
electronic state population of the atoms and molecules. The NEQAIR 3 code applies a 
collisional-radiative (CR), or non-Boltzmann, model for the calculation of the electronic state 
populations, while the RAD/EQUIL 4 code assumes a Boltzmann distribution of electronic states. 
This difference is significant because in regions of chemical and thermodynamic nonequilibrium, 
the electronic states of atoms and molecules do not follow a Boltzmann distribution. The non- 
Boltzmann model applied in NEQAIR is based on the approach originated by Bates et al. 5 . 

Figure 1 illustrates a typical population distribution of the electronic states for 
nonequilibrium atomic nitrogen. The populations were calculated using the CR model 
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developed in the present work and presented in Section VII. In this figure, V, is the number 
density of level i, gi is the degeneracy of the level, and E, is the term-energy of the level. It is 
seen that the CR model predicts much lower number densities for the highly excited states than 
predicted by the Boltzmann model. This is a characteristic of nonequilibrium conditions during 
compression, such as those found behind a shock wave in hypersonic flow. The opposite is true 
for expanding nonequilibrium conditions, such as those found in a favorable pressure-gradient 
boundary layer with chemically reacting flow. The lower excited state number densities result in 
decreased radiation, relative to the Boltzmann prediction, from flowfield regions of chemical 
nonequilibrium. 

In the early 1990s, flowfield calculations that included chemical and thermodynamic 
nonequilibrium became standard. For these flowfields, the RAD/EQUIL code was inadequate 
because it assumed a Boltzmann distribution of the electronic states, which is incorrect in regions 
of nonequilibrium. Instead of resorting to the computationally intensive NEQAIR code for these 
situations, corrective procedures and major revisions to the RAD/EQUIL code were introduced. 
Gaily et al. 6 ' 7 developed two approximate correction methods, which were implemented in 
RAD/EQUIL, to account for the non-Boltzmann population of atomic electronic states. The first 
of these models, referred to as the 1 st order local thermodynamic nonequilibrium (LTNE) model, 
was originally proposed by Carlson 8 . It assumed that the excited atomic states were in 
equilibrium with the ions and electrons, instead of with the ground electronic state of the atom. 
The excited state number densities were therefore calculated using the Saha-Boltzmann equation 
as follows 


N? b = N+N, 


l h 2 \ 


3/2 


2mnkT L , 

instead of with the regular Boltzmann equation 


cx p[- lu(E i -E ionize )/kT e ] 
2 < 2 + 


( 1 ) 
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( 2 ) 


Nf = N, 


g,exp(-hcE;/kT e ) 

Q a 


On the other hand, the number densities for the ground states, which include three low-lying 
states for nitrogen and oxygen, were calculated using Eq. (2). The rationale for such a simple 
model is apparent from Figure 1, where the results of Eqs. (1) and (2) are compared with the 
detailed CR model results. It is seen that the three low lying states do follow a Boltzmann 
distribution closely, and the highly excited states do approach the Saha-Boltzmann population as 
the ionization limit is approached. The 1 st order LTNE method is a very simple method for 
approximately accounting for the non-Boltzmann population of atomic states. It also shows the 
connection between chemical nonequilibrium and non-Boltzmann radiation. In regions of 
chemical equilibrium, Eqs. (1) and (2) are equivalent, by definition, because Eq. (1) is obtained 
by substituting Eq. (2) into the “law of mass action” equation 9 for the following reaction 


A A + +e~ (3) 

where A represents a neutral atom and A + represents its ion. In regions of chemical equilibrium, 
the flowfield solver will have enforced the law of mass action, which requires for the reaction in 
Eq. (3) that the neutral species, ion, and electron number densities are related as follows 


N„ = AEAf 


< h 2 N 
2mnkT„ 


3/2 


Q± 

2 £ 


-exp (hcE ionize lkT e ) 


(4) 


Substituting this into Eq. (2) makes clear the equality of Eqs. (1) and (2) in regions of chemical 
equilibrium and explains why, in such a case, all three of the lines in Figure 1 would be the same. 
A procedure similar to the 1 st order LTNE model was also outlined by Geendyke and Hartung 10 , 
although they applied the method exclusively to a step model representation of the spectrum. The 
2 nd order LTNE model, also presented by Gaily et al. 11 , grouped together the excited states of an 


atom and treated them as a single species in the flowfield calculation. The populations of the 
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various levels of the excited states were then distributed in a Boltzmann distribution, with the 
number density of the “excited” specie determined from the flowfield calculation. The 2 nd order 
model was shown to result in slightly larger values than the 1 st order model, with a difference of 
about 10% for the cases considered. 

III. Atomic Level Model 

The emission and absorption of an atomic line is proportional to the number density of the upper 
and lower electronic state of the defined transition, respectively. If the number density of these 
states are assumed to follow a Boltzmann distribution, and if the energy and degeneracy of the 
upper or lower states are defined for each atomic line, then there is no need for a separately 
defined set of atomic levels. This is true because the Boltzmann relationship and Plank function 
allow for the emission and absorption of each line to be determined. For the non-Boltzmann 
case, on the other hand, a separate model is required for calculating the population of the 
electronic states of the radiating atom. The definition of each atomic state in this calculation must 
be related to the defined upper and lower state of each atomic line transition. Therefore, in 
preparation for the discussion of the atomic non-Boltzmann model (in Section VII), this section 
presents a model for the atomic levels of nitrogen and oxygen suitable for the present hypersonic 
shock-layer applications. 

There are few available guidelines for defining an adequate set of atomic levels for 
nitrogen and oxygen for the purpose of modeling the radiative heat flux from a chemical 
nonequilibrium shock-layer. Bourdon and Vervisch ~ discussed the choice of an energy level 
model for atomic nitrogen. The focus of their study was the computation of the three-body 
recombination rate through the use of a CR model. They concluded that the 22 level model 
presented by Park 13 and the 14 level model presented by Kune and Soon 14 were inadequate and 
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proposed instead a 43 level model similar to that used by Park 15 . Bourdon et al. 16 reached a 
similar conclusion that 41 levels are required for atomic oxygen. Because the present study is 
focused on calculating the radiative flux from a shock layer, and not calculating the three-body 
recombination rate, the conclusion of Bourdon and Vervisch 12 and Bourdon et al. 16 that over 40 
levels are required for nitrogen and oxygen is not directly applicable. Because no further 
guidance is available, the current model will be chosen as a compromise between the Bourdon 
studies and the widely used (for shock layer radiation) Park 15 model. The new model, which is 
presented in detail by Johnston 18 , consists of 35 levels for N and 32 levels for O. The data for 
these levels were obtained from the NIST database 17 ; except for those that group all the levels of 
a single principal quantum number (n), which were taken from Park’s work. For both N and O, 
the first 25 levels are ungrouped, while the rest of the levels consist of groups of closely spaced 
levels. This allows for maximum precision in the CR model calculation for the lower levels, 
which are the most important for the radiation calculation. It is likely that the total number of 
levels could be reduced, with a negligible loss in accuracy, by further grouping. Nevertheless, 
the present groupings will be maintained to provide benchmark results for assessing further 
simplifications. The relationship between the energy for a grouped level i and a set of ungrouped 
levels i ’ is: 



and for the degeneracy is 

Si=^Sr ( 6 ) 
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where the summation is over all of the levels i’ in the group i. The levels contained in a group 
are assumed to be populated in a local Boltzmann distribution, which allows the individual level 
number density (Nr) to be related to the group number density (/V,-) as follows 


N r = Nf —exp 
8 : 


kT„ 


( 7 ) 


This relationship is useful for relating data given in terms of the number density of the individual 
level (such as atomic line data and electron-impact excitation rates) to the number density of the 
appropriate grouped level. 


IV. Excitation and De-Excitation Mechanisms 

As mentioned previously, the population of the radiating atomic and molecular electronic states 
required for the calculation of the radiative emission and absorption are calculated by solving the 
so-called Master Equation. The derivation of this equation requires the specification of the 
various repopulation and depopulation mechanisms of the atomic or molecular states. The 
present model will consider electron impact, heavy particle impact, and radiative processes. 
These processes are defined in Table 1 and discussed in detail by Johnston 18 . 


V. The Master Equation 

The Master Equation is the governing equation for the electronic state populations of atoms and 
molecules. This differential equation, which must be solved for every state of a radiating atom or 
molecule, equates the time-rate-of-change of a level’s population with all of the populating and 
depopulating mechanisms listed in Table 1. It is common to express the electronic state 
populations as a nondimensional population factor p„ which is defined as 


Pi = 


N; 


N, 


SB 


( 8 ) 
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SB 

where N,- is the number density of the electronic level i and N, is number density of level i 
predicted by the Saha-Boltzmann equation defined in Eq. (1). This factor is convenient because it 
forces the values from all levels of a given atom to have similar magnitudes. The Master 
Equation will be written in terms of this factor. Other non-dimensional quantities that will be 
useful in the present discussion are defined as follows 


iPi = 


K E 

n se 


(9) 


Xi = 



( 10 ) 


The Master Equation is derived by equating the time-rate-of-change of level i to the excitation 
and de-excitation processes presented in Table 1. This results in the following equation: 


dN t 

dt 




+ 

elec -imp -dis 



+ 

rad -ex 



rad -ion 


( 11 ) 


Eq. (11) represents a set of m differential equations for the m unknown p, values. If it is assumed 
that dNj / dt is equal to zero, then this set of differential equations reduces to a set of linear 


algebraic equations for the p, values. This is the quasi-steady state (QSS) assumption, which is 
known to be invalid for the ground state of atoms and molecules 19 . To accommodate this, the 
Master Equation may be solved using the QSS assumption for all levels except for the first. The 
population of the first level is then obtained by equating the sum of the level number densities of 
the species to the total number density of that species, as will be discussed in the next section. 

The solution of the Master Equation may be written in matrix form for all of the levels 
except for the first as follows 18 

Th.m = r 0 + r 0() ip + r,p, (12) 
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where 



K e (i,\)+ 08) 

With r Q0 = 0 , Eq. (12) is equivalent to Eq. (14) presented by Bates et al. 5 . Note that Eqs. (14) to 
(18) are functions of the following flowfield variables: Nx, Ny, Nm, N e , T, and T e (the T e 
dependence is implied through the rate coefficients for the electron-impact processes and T, for 
the heavy-particle processes). Although N+, N? B , and nP e appear in these equations, they may 
be reduced to a N e , N x , N Y , T t , and T e dependence using the Saha-Boltzmann equation (Eq. (1)) 
and detailed balancing 13 ' 18 . For atomic species, it is well known that the heavy-particle impact 
excitation does not contribute significantly to the excitation and de-excitation processes. Ignoring 
these terms in Eqs. (14) to (18) for atoms results in these equations being functions of only N e 
and T e (note that the dissociation processes are not considered for atoms). For molecular species, 
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it will be shown that for the lunar-return shock layers of present interest, the heavy particle 
excitation and electron-impact recombination processes may be ignored as well. This results in 
Eqs. (14) to (18) also being a function of only N e and T e . The significance of these dependencies 
is that they result in the r 0 and r. vectors being functions of only N e and T e , which means that 

these vectors may be conveniently compiled in tables or curve fits in terms of these two 
variables. Thus, the computationally intensive matrix inversion required to calculate these 
vectors, shown in Eq. (13), may be avoided. This was recognized by Bates and Kingston' , who 
present tables of r 0 and f\ for hydrogen plasmas with 0 < N e < oo ; 250 < T e < 64000, and A j t = 0 

(optically thick) and 1 (optically thin). The present curve-fits developed for the atoms and 
molecules of present interest will be discussed in Sections VII and VIII, respectively. 

VI. A Simple Method for Calculating p i 

The previous section presented a method of calculating the p, values for energy levels above the 
first level by making the quasi-steady state approximation ( dN i /dt= 0). Equation (12) shows 
that these values can be represented as a linear combination of pj. This section discusses a 
method of calculating pi using the following equation 

m 

2 n i-n. < l9 > 

7=1 

which forces the sum of the individual state populations to equal the total number density of the 
species obtained from the flowfield. This approach, which was introduced by Park' , is used in 
NEQAIR’s QSS subroutine 31922 . Rewriting Eq. (19) in terms of p, and separating pi from the 
summation results in 

m 

PiXi + 2 PjXj = 1 ( 2 °) 

J = 2. 
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Substituting Eq. (12) into this expression and solving for pi results in the following 


m 

1 _ 2 2 00+ F oo 



( 21 ) 


In summary, the solution procedure for this approach involves the following: 

1) Obtain the r 0 , r 00 , and r\ vectors either by solving Eq. (13) or from a previously prepared 
table or curve fit (recall that r 0 and f\ contain m terms). 

2) Calculate pj from Eq. (21) 

3) Calculate the p, values from Eq. (12), which is written for this case as 

Pi = r 0 (i )+ r m (i M0+ h if )p, (22) 

As mentioned previously, the method outlined here is equivalent to that used by Park in the QSS 
subroutine of NEQAIR, although in the present case the equations have been formulated in terms 
of the r Q , r 00 , and r x vectors. 


VII. Non-Boltzmann Modeling for N and O 

This section presents rate models, specifically for atomic nitrogen and oxygen, for the significant 
excitation processes presented in Table 1 and applied in Section V for the Master Equation. The 
various rates were collected from available experimental and theoretical published sources, 
compared, and a “best” value chosen. For some processes, no rates were available in the 
literature, and so approximate formulas were applied. The majority of the effort in constructing 
this model was obtaining electron-impact excitation rates for transitions between all of the levels 
of N and O. This effort will therefore receive the majority of the following discussion regarding 
the various excitation processes. 
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A) Electron-Impact Excitation: 

The process of electron impact excitation was discussed in general in Section IV. The fact that 
there are electron-impact excitation processes between each atomic level results in a large 
number of rates being required. For the present study, the 35 levels for N and 32 levels for O 
imply that there are 595 and 496 rates required for each species, respectively. Consequently, 
collecting this complete set of rates is a difficult and tedious task. The difficulty is compounded 
by the fact that there are a limited number of detailed calculations and experimental studies for N 
and O, which are furthermore limited to a relatively small number of transitions. Thus, extensive 
use of approximate analytic formulas is required to complete the set of rates for each species. 
These approximate formulas are discussed in the following paragraph, followed by a discussion 
of the available detailed rates for both N and O. 

A summary of approximate formulas for the electron-impact cross sections are presented 
in Table 2. Although these approximate formulas are relatively old, a review of recent 
astrophysics' ’ ’ ~ " and plasma physics ~ literature reveals that in the absence of 
experimental data or detailed calculations, these approximate formulas are the best available 
alternative. These approximate formulas will be compared with detailed calculations in the 
following paragraphs. This comparison is aimed at determining the best available rate for each 
transition, as well as confirming the validity of the approximate formulas for transitions with no 
detailed rates available. The detailed rate or approximate formula chosen for each transition is 
presented at the conclusion of this discussion. 

The two main sources of detailed quantum mechanical calculations for nitrogen are the 
recent works by Frost' and Tayal' . Frost' presents tabulated values for the transition rates from 
the lowest three levels (i = 1, 2, 3) to the first 21 levels. These rates were obtained using an R- 
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matrix approach and were shown to match available experimental data. Tayal' presents 
tabulated values for the transition rates from the lowest three levels (i = 1, 2, 3) to the first 12 
levels. An R-matrix approach was also used in this study, although comparisons with 
experimental results or other calculations were not made. Figure 2 compares the calculations of 
Frost and Tayal with the approximate formulas summarized in Table 2. Figure 2(a) shows that 
for the 1-5 forbidden transition, the results of Frost and Tayal differ by an order of magnitude, 
while the approximate expressions are all similar and fall between Frost’s and Tayal’s results. 
Figure 2(b) shows that for transitions between higher levels (3-20), Allen’s approximate 
formula provides values nearest the values predicted by Frost. 

From the above comparison, along with the study of many transitions not discussed 
above, it was concluded that Van Regmorter’s expression is the best available approximate 
expression for allowed transitions where the upper level of the transition is less than j = 22. 
Similarly, Allen’s formula was concluded to be the best for forbidden transitions with j less than 
22. For transitions with j greater than 22, Gryzinski’s formula is accurate because of the 
hydrogenic nature of these upper levels, and was therefore applied. A summary of these chosen 
rate sources for the final nitrogen model is presented in Table 3. Note that Frost’s detailed 
calculation values were applied where available. These were chosen instead of Tayal’s values 
because Frost considered more transitions, including all of those considered by Tayal. 

The main sources of detailed electron-impact excitation rates for atomic oxygen are the 
studies by Laher and Gilmore , Bhatia and Kastner , and Zatsarinny and Tayal . Laher and 
Gilmore ' obtained rates for transitions from the ground state to the first 19 excited states by 
reviewing past experimental data and by using a correlation presented by Jackman et al. 34 . Bhatia 
and Kastner' used a distorted-wave calculation to obtain rates for most of the transitions 
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between the lowest ten states. Some of the rates presented by Bhatia and Kastner were taken 
from Laher and Gilmore, although they were not specified. Zatsarinny and Tayal calculated 
values for transitions from the lowest 3 states to the first 21 excited states using an R- matrix 
approach (with many values missing for transitions from the second and third state). Table 4 lists 
the sources of the rates chosen for the final oxygen model. 

B) Electron-Impact Ionization: 

The process of electron impact ionization is of less importance than electron-impact excitation 
and bound-bound radiative transitions, and therefore it does not need to be treated as precisely. 
Also, only one rate is required for each level so that significantly fewer rates are required than 
the bound-bound processes. 

The rate coefficients for electron-impact ionization from excited levels may be calculated 

o c 

accurately with the following formula proposed by Drawin 

2 

^(/,c)=1.46xlO“ 10 r e 1/2 ^ fyV'iCy) (roV 1 ) (23) 

v ^ ionize , 

where § is the number of equivalent electrons ( § = 3 for the ground state of atomic nitrogen and 
oxygen and § = 1 for all excited states) and y is the reduced energy of the incoming 
electrons, y = {E ionize - E i )hc/kT e . The function ip] is approximated by the following equation 

T7T— + ln [l.25(l + -)|} (24) 

1 + y [ 20 +y [ ( y jjj 

For ionization from the lowest two states of oxygen and nitrogen, the rate coefficients proposed 
by Kune and Soon 14 and Soon and Kune 36 were used. 

C) Bound-Bound Radiative Transitions (Atomic Line Transitions): 
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The influence of bound-bound radiative transitions on the Master Equation was mentioned in 


Section IV. This process was also discussed in Section II of Ref. 1, as it related to the spectral 
modeling of atomic line radiation. The set of atomic lines for nitrogen and oxygen chosen from 
the NIST database in Ref. 1 are also applied in the present Master Equation model. In general, 
the transition probabilities are related to the oscillator strengths as follows 




%n 2 e' > 




me 


i,j °j 


gf ( 

he ( 

17 ^ 

— -exp 


'/ - Ej J 

8} l 

kT e V 



(V 1 ) (25) 


where the term in brackets accounts for grouped levels, as discussed in Section III. The 
summation in this equation is over all of the i’ and /’ individual levels present in the i and j 
grouped levels. Only optically allowed radiative transitions are considered in the present study. 
Although Kune and Soon 14 and Soon and Kune 36 considered various forbidden transitions, they 
note that their influence is small, especially the relatively large electron number densities present 
in the hypersonic shock-layers of interest in this study. 

D) Bound-Free Radiative Transitions: 

The bound-free radiative transitions, introduced in Section IV, are of minor importance for most 
cases. The approximate expression presented by Drawin , based on the hydrogenic model for 
the bound-free cross sections, is applied in this study. This is written as follows: 

A(c,/)= 5. 20xl0“ 14 Zw 1 3/2 G ! r 3 .£■(«, )exp(w ; ) ipn 3 s _1 ) (26) 


where the exponential integral is defined as 

*(„,)- j exp H , (27) 

t=Uj 

The detailed bound- free cross-section obtained from the TOPbase could be applied for a more 
detailed calculation, although this has a negligible influence on the results. 
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VIII. The Approximate Atomic CR (AACR) Model 


Using the detailed CR model described in the previous section as a baseline, an approximate 
atomic CR (AACR) model was developed, which is significantly simpler to apply than the 
detailed CR model and is nearly as accurate at predicting the resulting nonequilibrium radiative 
emission. The details of this method are presented here and its results compared with the detailed 
CR model for some relevant shock-layer conditions. Comparisons will also be made to Gaily ’s 
1 order LTNE method , which is the approximate method discussed previously in Section II. 
Recall that the Gaily’ s method assumed that the three lowest levels were in a Boltzmann 
distribution with the ground state (Eq. (2)) and the upper levels were in a Saha-Boltzmann 
distribution (Eq. (1)). 

The approximate atomic CR (AACR) model developed in this study contains the 
following three main approximations: 1) the r Q and f[ vectors defined in Eq. (13) are modeled as 
single curve fits over a wide range of T e and N e values, 2) closely spaced atomic states are 
assumed to be in a Boltzmann distribution with each other, 3) the three lowest atomic states are 
assumed populated by a Boltzmann distribution. The only approximate aspect of the first of these 
approximations is that a single curve fit over a range of T e and N e values deviates slightly from 
the actual curve-fitted data. The dependence of r Q and r x on only T e and N e was pointed out in 
Section V, and is therefore not an approximation. The second approximation listed above is 
apparent from Figure 1, which shows groups of closely spaced levels with nearly the same slope 
(among the levels in the group) as the Boltzmann and Saha-Boltzmann curves. With this 
assumption, once the population of a single level in each of these groups is known, then the 
population of the other levels in the group may be calculated from Eq. (7). Table 5 presents the 
grouping of levels for N and O, defined in terms of the level model discussed in Section III, that 
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were chosen here for the AACR model. The level in each group whose r 0 and r\ values were 
curve-fit are also listed in this table. These curve-fits were created from the detailed model 
discussed the previous section for the following range of conditions: 7,000< T e <14,000 K and 
lxl0 14 < N e <lxl0 16 particles/cm 3 . The third assumption that the first three levels are populated 
by a Boltzmann distribution, which can be shown to be a good approximation for all the 
conditions of interest in the present study, allows the number densities of the curve-fit level for 
each group to be written in terms of Eqs. (8) and (22). The procedure for applying the AACR 
method is as follows: 

1) Calculate the r 0 and r/ values for each of the groups listed in Table 5 (except the first 
group) from the curve fits presented in Appendix C of Ref. 18. These ro and n values 
belong specifically to the levels listed in the column “Curve-Fit Level for the Group... ” . 

2) Calculate the number densities for the levels in the first group listed in Table 5 assuming 
a Boltzmann distribution (Eq. (2)). 

3) Calculate the number densities of each curve-fit level using the ro and rj values calculated 
in Step 1, the ground state number density from Step 2, and Eqs. (8) and (22). 

4) Calculate the number densities for the other levels in each group using Eq. (40) and the 
number densities calculated in Step 3. 

To examine the accuracy of the proposed AACR model relative to the detailed atomic CR model 
and the approximate model proposed by Gaily, the radiation for the Fire II 1634 second was 
calculated using each of these models, along with the radiation properties listed in Ref. 18 
(which are slightly different than those in Ref. 1) and an uncoupled LAURA stagnation-line 
flowfield. The temperature profile for this flowfield is presented in Figure 3. This figure 
indicates the large region of nonequilibrium present in the shock-layer. The total wall-directed 
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radiative flux values resulting from each of the atomic non-Boltzmann models are presented in 
Figure 4. Also presented is a case where the atoms were assumed to be populated in a Boltzmann 
distribution. The molecular species are treated in each case using the non-Boltzmann model 
presented later in this paper. Figure 4 illustrates the significant result that three CR models 
predict very similar radiative flux values throughout the shock layer, with the AACR and Gaily 
Model within 2% and 5% of the detailed CR model, respectively. Although the AACR model 
produces better agreement with the detailed CR model than the Gaily model, it should be noted 
that the Gaily model is significantly simpler conceptually and in application. The importance of 
applying a non-Boltzmann model, rather than assuming a Boltzmann distribution, is clearly 
shown in this figure with the nearly 100% over prediction of the Boltzmann case relative to the 
various CR model cases. The excessive radiation from the Boltzmann model originates entirely 
from the region thermochemical nonequilibrium directly behind the shock. The three CR model 
cases, on the other hand, show that this region contributes only slightly to the radiative flux. 

The good agreement between the radiative flux predicted by the two approximate models 
with that predicted by the detailed CR model was noted in the previous paragraph for the Fire 
1634 s case. To investigate this comparison further, two points in the shock layer were studied in 
detail: a nearly thermochemical equilibrium point at z = 2 cm and a nonequilibrium point near 
the shock at z = 4 cm. The temperature and relevant number densities for these points are 
presented in Table 6. 

For the nonequilibrium point at z = 4 cm, Figure 5(a) compares the number densities, 
divided by the degeneracy, of the excited states predicted by the three models (the lowest three 
levels are not shown because the result of all three models are indistinguishable and the values 
are orders of magnitude larger than those shown). The significant under prediction of the excited 
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state number densities by the Gaily model is apparent in this figure, as well as the excellent 
agreement of the AACR model with the detailed CR model. The under prediction by Gaily ’s 
model could have been predicted from Figure 1, which presents the excited state number 
densities for similar conditions and illustrates the trend that the Saha-Boltzmann values for the 
excited states, which in Gaily’ s model are assumed equal to the excited state number densities, 
are lower than the CR model prediction. Figure 5(b) presents the nitrogen atomic line emission 
from each level corresponding to the number densities shown in Figure 5(a). The differences in 
these emission values are consistent with the differences in the number densities of these levels, 
as they should be since the emission is proportional to the number density of the upper level of 
the transition. For the detailed CR model case, the total emission from nitrogen lines is 0.065 
W/cnr-sr while the overall emission from all radiative processes is 0.844 W/cnr-sr. Thus, the 
nitrogen lines contribute less than 10% to the emission (most of the emission at this point is due 
to molecular band systems), which explains why the radiative flux at z = 4 cm in Figure 4 
predicted by the Gaily model is similar to the other models, even though it was shown in Figure 
5(b) that this model significantly under predicted the nitrogen line emission. 

The z = 2 cm point listed in Table 6 is essentially in chemical equilibrium. This means 
that the Gaily model reduces to a Boltzmann distribution (because the Saha-Boltzmann and 
Boltzmann equations are identical for chemical equilibrium conditions). Figure 5(c) compares 
the excited states number densities for this case, and it is seen that the population of the lower 
levels predicted by the Detailed and AACR model are noticeably lower than the Gaily 
prediction. These levels do not reduce to their Boltzmann values because they are strongly 
radiating, meaning they are rapidly depopulated by spontaneous emission of photons and 
corresponding transitions to lower levels. The characteristic of being a strongly radiating level is 


22 

American Institute of Aeronautics and Astronautics 



a result of being the upper level of atomic line transitions with large transition probabilities, 
which were defined in Eq. (25). Also, it has been assumed that the gas is optically thin, as it 
applies to the Master Equation, meaning the escape factors have been set equal to one. This 
maximizes the influence of radiative transitions on the Master Equation. Figure 5(d) compares 
the nitrogen atomic line radiation for this case. As expected from the number densities, the Gaily 
model over predicts the emission for the lower states, which are seen to be the most strongly 
radiating. The total (including all radiating mechanisms of all species) emission from this point 
for the Detailed CR model case is 2.29 W/cnr-sr while that from nitrogen line emission is 1.02 
W/cnr-sr. Thus, it is seen that the nitrogen lines contribute a significant fraction of the total 
emission at this point, unlike for the previous z = 4 cm point. The Gaily model predicts a 1.22 
W/cnr-sr contribution from the nitrogen lines. This larger value explains the larger radiative flux 
value from the Gaily model shown in Figure 4. 

From the present discussion regarding the non-Boltzmann modeling of atomic species, it 
may be concluded that for the shock layer applications of present interest, the AACR model 
provides a sufficiently accurate non-Boltzmann model that is both computationally efficient and 
simple to apply. The six grouped levels used by this model for each species (listed in Table 5) 
essentially means that the 35 and 32 levels used to model N and O may be reduced to these six 
grouped levels. 

IX. Non-Boltzmann Modeling for N2 + 

For shock layers in air at lunar return conditions and near peak heating (radiative and 
convective), the N 2 + and N 2 molecules are the main contributors to the nonequilibrium molecular 
band radiation. The contributions from the NO and 0 2 bands discussed in Section VI of Ref. 1 
are limited to the equilibrium regions of the shock-layer. For lower velocity nonequilibrium 
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conditions, such as those studied by Levin et al. , the shock layer temperature is on the order of 
6,000 K and the nonequilibrium contribution from NO and O 2 may be large. For such conditions, 
the present discussion is not applicable, and instead, the discussions by Levin et al. 37 and 
Gorelov et al. should be consulted. The present section discusses the modeling of only N? , 
while a similar discussion is presented by Johnston 18 for N 2 . 

The rate model for the electronic excitation of N 2 + , for the lunar return shock layers of 
interest here, does not require that all of the processes discussed in Section IV be included. 
Nevertheless, for completeness, all of the processes will be included in the initial model 
proposed here, except for the electron-impact ionization process, which was shown to be 
negligible in the preliminary stage of this study (using the rates presented by Teulet 39 ). Note that 
the process of electron-impact ionization of molecules is not considered in the kinetic model 
applied for the flowfield chemistry 13 . 

The electronic levels for N 2 + are shown in Table 7 and the chosen rates for electron- 
impact excitation, dissociation, heavy-particle excitation, and radiative transitions are shown in 
Table 8 along with their reference source. The rates are related to the coefficients listed in this 
table as follows: 

K r = ct r T e br exp(- c r /T e ) (28) 

where K, represents K e (iJ ) for r = 1 to 6, K e (i,d) for r = 7 to 10, and K M (iJ ) for r = 11 and 12. For 

IQ 

many of the rates, the only available values were the theoretical predictions by Teulet et al. . 
These values were used when no other values were available. Preference was given to other 
values because the method used by Teulet et al. was more approximate than the experimental and 
detailed calculation procedures used by other researchers. The most important rate, and 
fortunately the most frequently studied rate, is that for the electron impact excitation from the 
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N 2 + (X) state to the N 2 + (B) state, which is the second rate listed in Table 8. This rate is the most 
significant for determining the population of the N 2 + (B) state, which is the upper level of the 
important N 2 + (l-) band system. Figure 6 compares the rates proposed by numerous researchers 
for this process. The solid lines represent the cases where the actual rates were presented by the 
researcher, whereas the non-solid lines indicate that the excitation cross-section was presented, 
which was converted to a rate by numerically integrating the following formula: 



This figure indicates the large discrepancy between the various studies, with values ranging 
many orders of magnitude. Because there is no guidance in choosing the best of these values, the 
fixed-nuclei result of Nagy 40 was chosen because it was the most recent and because it was 
roughly in the middle of the other predictions. The other sources of rates used for N 2 + in the 
present study are Gorelov et al. , Flagan and Appleton 41 , and Chernyi and Losev 42 . 

The first topic to examine for the proposed N 2 + model is the importance of each of the 
processes listed in Table 8. To aid in this study, it is useful to define the ratio of the inflow and 
outflow of a level j due to a specific process. For example, the ratio of the outflow from level j 
due to electron-impact transitions to the total outflow from level j (due to all process) may be 
written as follows: 

m 

J / K e (j,i)v j N e 

{ratio of j outflow due to el - imp - trans } = — ( 30 ) 

\ ^ J outflow 

where the total outflow from level j is written as 
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The equations for the outflow and inflow from each of the processes discussed in Section IV may 
be written analogously. We are interested in the emission from the N 2 + (l-) band system, which 
emits from the N 2 + (B) state, which for the present purposes represents j in Eqs. (30) and (31). 
Figures 7 and 8 present the inflow and outflow, respectively, for the Fire II 1634 s case using the 
rates listed in Table 8. The relative contributors to the inflow, or transition of molecules into the 
N 2 + (B) state, are shown in Figure 7 to be dominated by electron-impact transitions from the 
N 2 + (X) state, which is the important process discussed previously (r = 2 in Table 8). The other 
electron-impact transition processes are seen to be much smaller, although not negligible 
throughout the shock layer. The electron-impact dissociation and heavy-particle impact 
transitions are seen to only contribute noticeably in narrow regions near the shock and the body. 
The same trend is apparent in Figure 8, which shows the relative contributors to the outflow, or 
transition of molecules out of the N 2 + (B) state. The influence of radiative transitions is seen to 
contribute significantly to the outflow. Note that this calculation assumes that the molecular 
radiation is optically thin, which means that the radiation transitions have their maximum 
possible influence on the outflow. Figure 9 presents the net inflow, meaning the inflow ratio of 
Figure 7 minus the outflow ratio of Figure 8, for the present case. Because the quasi-steady state 
assumption was applied is solving the Master Equation, meaning the right-hand side of Eq. (11) 
was set to zero, the sum of the various net inflow components equals zero at each point in the 
shock layer. This figure shows that although the electron-impact inflow and outflow transitions 
to and from N 2 + (A) and N 2 + (C) have noticeable components in Figures 7 and 8, they are of equal 
magnitude and cancel out to a negligible amount for the net magnitude shown in Figure 9. 
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Therefore, they have little influence on the solution of the Master Equation and the resulting state 
populations. 

The conclusion that the electron-impact dissociation and heavy-particle impact transitions 
do not contribute significantly to the population of the N 2 + (B) state is applicable to all of the 
shock layers examined in this study. If these mechanisms are ignored, then from Eq. (12), the 
solution of the Master Equation becomes dependent upon only the electron number density and 
the electronic temperature, as mentioned previously. This dependence on only these two 
variables allows for a curve-fit to be conveniently constructed for the Ni + state populations. 
These curve-fits were constructed using the present rate model and are presented in Appendix C 
of Ref. 18 for the four N 2 + levels listed in Table 7. Unlike the curve-fits for the atomic species, 
which consisted of fits for ro and r/, the number densities of the levels are fit directly for 
molecules. This is possible for molecules because r 0 is equal to zero, which as seen from Eqs. 
(16) and (13) is a result of the electron-impact ionization and bound- free radiation being ignored. 

A comparison of the number densities for N 2 + obtained with the CR model presented here 
and those obtained with the Boltzmann model are presented in Figure 10 for the Fire II 1634 s 
case. The largest deviation from a Boltzmann distribution is seen to be present for the N 2 + (B) 
state. The CR model predicts a much lower number density, which implies that the emission 
from the N 2 + (l-) band is reduced proportionately. Figure 11 illustrated this reduction in emission 
by comparing the wall directed radiative flux profiles (resulting from all radiation mechanisms) 
for the Fire II 1634 s case obtained by assuming a Boltzmann distribution of the N 2 + (B) state or 
by using the present CR model for this state. This figure also shows the effect of choosing other 
rate values for the r = 2 reaction listed in Figure 6. It is seen that the Teulet et al. value, which 
was the largest of the rates considered, drives the N 2 + (B) state to a Boltzmann distribution, while 
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the Gorelov et al. value, which was the smallest of the rates considered, decreases the flux 
about 10% below the present rate value. 

X. Conclusions 

A simplified method of solving the Master Equation for atomic and molecular electronic state 
populations was presented. A novel approach of curve-fitting the non-Boltzmann population of 
the radiating atomic and molecular states as a function of the electronic temperature and electron 
number density was developed. The approximate atomic CR model (AACR) was developed 
based on this curve-fit treatment. This approach, which reduced the number of atomic levels 
treated, was shown to provide electronic state number densities within 5% of the detailed model. 
A new set of excitation rates was compiled for the significant radiating species in a lunar-return 
shock-layer, which include N, O, and N 2 + . The influence of electron-impact dissociation and 
heavy-particle impact transitions were found to be negligible for the N 2 + molecule. The 
difference between the rate values obtained from various sources was shown, and influence of 
these values on the resulting radiative flux for the Fire 1634 s case was illustrated. 
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Tables 


Table L Excitation processes accounted for in the present treatment of the Master Equation 
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Table 2. Summary of Approximate Equations for Electron-Impact Excitation Rates 



Name 

Transition Type 

Eq. in Ref. 18 

Ref. 

1 

Van Regmorter 

allowed 

(4.30) 

43 

2 

Park approx. 

forbidden 

(4.30) times 0.4 

44 

3 

Gryzinski 

any 

(4.32) -(4.34) 

31 

4 

Mien 

forbidden 

(4.31) 

30 


Table 3. Selection of Electron-Impact Excitation Rates for N 

1) Frost’s rates are used for i =1,2,3 and j = 2 through 21 

2) Remaining allowed transitions with j < 22 use Van Regmorter’s formula 

3) Rem ain i ngforhidden transitions with j < 22 use Allen’s formula 

4) All remaining transitions use Gryzinski’s formula 


Table 4. Selection of Electron-Impact Excitation Rates for O 

1) Zatsarinny and Tayal rates are used for i =1,2,3 and j = 2 through 21 where available 

2) Bhatia and Kastner for remaining rates with j < 10 where available 

3) Gordillo and Kune for remaining rates withy < 7 

4) Remaining allowed transitions withy < 22 use Van Regmorter’s formula where available 

5) Remaining forbidden transitions withy < 22 use Allen’s formula 

6) All remaining transitions use Gryzinski’s formula 



Table 5. Level Groupings for the AACR Model 


Group # 

Range of Levels for 

Curve-Fit 

Range of Levels for 

Curve-Fit 


N from Table 3.1 of 

Level for the 

O from Table 3.2 

Level for the 


Ref 18 

Group for N 

of Ref. 18 

Group for O 

1 

1-3 

- 

1-3 

- 

2 

4-6 

5 

4-7 

5 

3 

7-13 

10 

8-13 

10 

4 

14-21 

17 

14-21 

17 

5 

22-27 

24 

22-27 

24 

6 

28-35 

28 

28-35 

28 


Table 6. Conditions for two points in the Fire 1634 s shock-layer 



z = 2 cm 

z = 4 cm 



Tel K) 

10,467 

12,670 



T'( K) 

10,459 

19,186 



2 

number densities (#/cm ) 




N 

1.74xl0 16 

1.2xl0 16 



N + 

3.90xl0 15 

3.5xl0 14 



e 

4.55xl0 15 

4.7xl0 14 



Table 7. Nj" 1 " electronic levels treated in the present model 39 

Level E e (cm 1 ) g e 
X 2 X g + 0.0 2 

A 2 n u 9,016.4 4 

B 2 Z u + 25,461.5 2 

C 2 I U + 64,542.0 2 



Table 8. Non-Boltzmann rate model for N 2 + 


r 

Process 

Cl r 

b r 

c r 

Ref. 

i 

N 2 + (x)+e" «> A(A)+e“ 

7.1xl0~ n 

0.0 

13,300 

38 

2 

N*_{x)+e~ jV 2 + (fi)+e“ 

2.0xl0~ n 

0.73 

36,649 

40 

3 

N*_{x)+e~ ^N 2 (c)+e" 

6.6xl0' 9 

0.41 

85,038 

39 

4 

AN(a)+<T ^N*{B)+e- 

1.0x10 9 

0.0 

23,500 

38 

5 

N 2 (A)+e“ ^N*{c)+e~ 

1.3xl0~ 7 

0.11 

78,403 

39 

6 

AN (fi)+<r <^N 2 (c)+e“ 

3.9x10 9 

0.34 

49,622 

39 

7 

N 2 (x)+<r ^ N + A +<T 

8.02x10 71 

5.54 

101,117 

39 

8 

N 2 (A)+e" o N + A +<r 

8.27x10 26 

4.38 

88,142 

39 

9 

N 2 (B)+e" N + A +e" 

2.58xl0 32 

5.81 

64,328 

39 

10 

N 2 (c)+e" «> N + A +<?" 

1.31 xlO 28 

4.93 

35,906 

39 


11 

n 2 (x )+ x 2 ^ n 2 (a)+ x 2 

3.8 xlO' 2 

-2.33 

12,978 

40 

12 

x 2 (x )+ N 2 o n* ( b)+n 2 

1.9 xlO' 2 

-2.33 

36,600 

41 

A- (s' 1 ) 

13 

N 2 (a)^ N 2 (x)+/jv 

6.7xl0 4 



42 

14 

N 2 (b)^ N 2 (x)+/!V 

1.5xl0 7 



38 

15 

n 2 (c)^ N 2 (x)+/iv 

1.4xl0 7 
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Figures 



Figure 1 . Population of the electronic states of atomic nitrogen in nonequilibrium 
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a) Forbidden transition between low N levels 
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b) Forbidden transition between a high and low N level 


Figure 2. Comparison of the various approximate electron-impact excitation cross-sections with 

detailed calculations. 



Figure 3. Temperature profile along the stagnation-line predicted by the LAURA code for the 

Fire 1634 s case. 



Figure 4. Influence of the various atomic CR models on the total wall-directed radiative flux for 
the Fire 1634 s case, the molecules are treated in all cases with the non-Boltzmann model 

presented later in this paper. 
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a) Number densities at z = 4 cm 



b) Line emission at z = 4 cm 
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c) Number densities at z = 2 cm 



Level Index, i 

d) Line emission at z = 2 cm 

Figure 5. Number densities and line emission from atomic nitrogen levels predicted by the 

various models. 
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Figure 6. Comparison of electron-impact excitation rates for populating the N 2 + (B) state. 
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Figure 7. The ratio of N 2 + ( B ) inflow for various process 
for the Fire II 1634 s case, using the Table 8 rates. 
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Figure 8. The ratio of N 2 + (B) outflow for various process 
for the Fire II 1634 s case, using the Table 8 rates. 



Figure 9. The ratio of N 2 + (B) inflow ratio minus the outflow ratio for each process corresponding 

to Figures 7 and 8. 
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Figure 10. Comparison of the electronic state number densities predicted by the 
Boltzmann and CR Model along the stagnation line for the Fire II 1634 s case 



Figure 11. Wall directed radiative flux profiles for the Fire II 1634 case 
for different treatments of the N 2 + (B) state population 


